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Abstract 



We study the retarded response function in scalar c/> 4 -theory at finite temperature. 
We find that in the high-temperature limit the imaginary part of the self-energy 
is given by the classical theory to leading order in the coupling. In particular the 
plasmon damping rate is a purely classical effect to leading order, as shown by 
Aarts and Smit. The dominant contribution to Landau damping is given by the 
propagation of classical fields in a heat bath of non-interacting fields. 



In the high-temperature phase of the standard model of electroweak interactions 
baryon-number and lepton-number violating processes are in thermal equilibrium [I]. 
This is of crucial importance for the presently observed cosmological baryon asymmetry 
and the 'sphaleron rate' of baryon-number violation is a corner stone of the theory of 
baryogenesis. 

In order to obtain the sphaleron rate one has to compute a real-time correlation 
function at finite temperature in a non-abelian gauge theory - a problem which, up 
to now, could not be solved. Several years ago Grigoriev and Rubakov suggested that 
the dominant contribution to the sphaleron rate could be obtained by computing the 
classical time evolution of the gauge fields and averaging over the initial conditions with 
the Boltzmann weight factor |J. This classical real-time method has indeed led to an 
accurate determination of the 'classical' sphaleron rate which was shown to be insensitive 
to the ultraviolet cutoff ||. The same procedure has been used to study other real-time 
properties of the SU(2)-Higgs model at finite temperature 

The classical real-time method has been questioned for several reasons. One worry 
concerns the well know ultraviolet divergencies of classical statistical field theory || [7| . 
Another open problem is the role of Landau damping §, ||, i.e., the loss of energy of 
quasi-particle excitations in the plasma to the heat bath. 

To elucidate the role of classical field theory for finite-temperature quantum field 
theory, Aarts and Smit have recently studied the plasmon damping rate in scalar 4 - 
theory ||10|| . As they have shown, the damping rate is determined by the classical theory 



to leading order in the coupling. In this paper we extend this analysis by studying the 
retarded response function in scalar 4 -theory. It turns out that the imaginary part of 
the self-energy is entirely given by the classical theory in the high-temperature limit to 
leading order in the coupling. We also calculate the dominant subleading contribution 
to the plasmon damping rate which is a genuine quantum effect. 

Retarded Response Function in the Quantum Theory 

We consider scalar 4 -theory at temperature T = 1/(3. Perturbation theory is most 
easily carried out in the imaginary time formalism which is based on the action 

S = J dx {^d^d^ + ^m 2 2 + ^A0 4 + 5) , 



dx = [ dr I d 3 x. (1) 

1/3 JO J 

Here m is the plasmon mass, 

- 2 = ^AT 2 , (2) 
and S stands for counter terms associated with zero-temperature ultraviolet divergencies 



as well as the resummation at finite temperature [ED]. The two-point function 



D(n - r 2 , xi - x 2 ) = (0(ti, xi)0(r 2 , x 2 )) (3) 



2 



satisfies the periodic boundary condition D(0, x) = D(j3, x). Hence, its Fourier transform 
is only defined for the discrete energies given by the Matsubara frequencies oj n = 27m//?, 
with n = 0, ±1, ±2, . . ., 



D(t,x) 



1 



d z x^D{iu n ,p). 



(4) 



Analytic continuation yields the function D(z,p) for complex values of z |14|]. This func- 
tion is related to the self-energy by a Dyson-Schwinger equation, 



Tl(z,p) = D- 1 (z,p)-D 1 (z,p). 



(5) 



Here D (z,p) = 1/ (z 2 — u 2 ), with uo p = y/p 2 + m 2 , is the free two-point function. Along 
the real axis real and imaginary part of the self-energy are defined by (u = ou*), 



U(u ±ie,p) = U R (uj,p) ± i T(u,p) . 



(6) 



Knowing D(z,p) in the complex plane, the various two-point functions can be eval- 
uated by choosing the appropriate path. We are particularly interested in the retarded 
two-point function which is given by 

/+oo A, i 
^ ^-D(u + i e,p)e~^. (7) 

For small couplings A the self-energy U(z,p) can be evaluated in perturbation theory. If 
D(z,p) has a pole at 

z = u> p -i%, (8) 
with | 

Q p- \J u l- n «K'P) > 7 P -^ r K ; p) ; (9) 

the integral in Eq. (|^) is easily evaluated, and one obtains the ordinary retarded Green's 
function, 

D%{t,p) ~ -6(t)— sin(u; p t)e-^* . (10) 

The retarded two-point function can also be expressed as ensemble average of a 
product of field operators, 

iD%(t -t',x- x' ) = Tr (e- pA [i(t, x), x 1 )]) Q(t - , (11) 

where H is the hamilton operator of the theory. Hence, Dr describes the response of the 
system to a perturbation caused by an external current. Adding to the hamiltonian a 
source term, 

H s (t) = J d 3 xj(t,x)4>(t,x), (12) 
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leads to a non-vanishing ensemble average of the field operator. In linear response theory 
one obtains 

($(t,x)) qu = Tr (e-^(t, £)) 

dt' I d 3 x'j(t', x')D R (t -t',x- x') . (13) 



The finite-temperature two-point functions have been studied in detail in the liter- 
ature. The self energy Il(z,p) has been calculated in perturbation theory to two-loop 



order fL2| , |13|j . We will be particularly interested in the damping rate, i.e., the imaginary 
part. It consists of two parts flUT] : the first one, which is present also at zero temperature, 
is due to the off-shell decay of a scalar particle into three on-shell scalar particles; the 
second part results from Landau damping, i.e., the scattering of the off-shell scalar par- 
ticle from the heat bath. The imaginary part of the self-energy reads at two-loop order 



(u > o) in, 

T(lu,p) = ^\ 2 (e^-l)jd$®[8{u-u> x -u, 2 -ui)hhh 

+ 35(oo + Ul -u 2 - w 3 )(l + /O/2/3] , (14) 
T(-oj,p) = -I>,p), (15) 

where d$>(p) is the integration measure, 

*® = (^(^(2^< 2 *> 3 ^- * " * " « • < 16 > 

and / is the Bose-Einstein distribution function, 

/H = ^n- ( 17 ) 

The two contributions to V(u, 0) have thresholds at uo = 3m and u = m, respectively. 
Only the first term is present at zero temperature. At temperatures T ^> m both terms 
are dominated by interactions with the heat bath, i.e., by Landau damping. 

Retarded Response Function in the Classical Theory 



The retarded response function can also be studied in the classical theory at fi- 
nite temperature. The classical hamiltonian is the sum of the free hamiltonian, a self- 
interaction term and a source term, 

H = Jd s x (V + l -(d<pf + \m\^ + iA0 4 + jcP) , (18) 

where j(t, x) is an external souce. Note, that m c i = l/l c i is an inverse length. It is different 
from the plasmon mass m given in Eq. (0) . The ensemble average of the classical field is 
then given by 

(</>(*, x)) d = || P7rP0e- H ^)0(t, x) , (19) 
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with 

Z = I VTrV^e' 11 ^ . (20) 
Here <p(t, x) is the solution of the equations of motion 

4>{t,x) = 7r(t,x), 

n{t,x) = A0(t,f)-m^(t,f)-^ 3 (t,f)-i(t,f), (21) 
with the initial conditions 

<f)(t , x) = 0(f) , ir(t , x) = ir(x) . (22) 

In Eqs. (0) and (|20"D the domain of integration is the space of initial conditions. 

The solution of the equations of motion satisfies the integral equation (t, t' > t , x = 

(*,£)), 

4>{x-j) = <j> {x) + J d 4 x'D R (x - x') j) + j(x')) . (23) 

Here <fio is the solution of the free equations of motion which satisfies the boundary 
condition, and Dr is the retarded Green's function, 

D R (t, x) = -J ^^e(t)i- sinKt) . (24) 

The classical solution can be expanded in powers of the external current. The term 
linear in j(x) is given by (cf. ([IB])) 



cj>{x) = / d A x'j(x')H R (x-x'), (25) 

where 



H R (x-x>)- 6 ^ 3: 



(26) 

j=0 



5j{x>) 

depends on the initial conditions <f)(x) and tt(x). From Eqs. (|23D and (|26"D one reads off 
the integral equation which determines H R , 

H R (x - x') = D R (x - x') d 4 yD R (x - y)<p 2 (y; 0)H R (y - x') . (27) 

In order to obtain the finite-temperature retarded response function, i.e., the classical 
analogue of D q R (cf. (|TT1)), one has to compute the ensemble average with respect to the 
initial conditions. This yields 



{4>{x)) d = J d 4 x'j(x')Di(x-x'), (28) 

where 

Dl\(x - x 1 ) = (H R (x - x')) d = ^J V(j)V7re^ H H R (x - x') . (29) 



Figure 1: The response function Hr expanded in powers of the classical field. Full lines 
denote retarded Green's functions Gr 



The classical retarded response function Dr can be evaluated by first expanding the 
solution (p7|) for Hr in powers of A (cf. Fig. 1) and then performing the thermal average 
for each term. This requires the evaluation of the thermal n-point functions 

(<p 2 (x 1 )...(P 2 (x 2n )) d . (30) 

Consider first the thermal average with the free Hamiltonian Hq. The corresponding 



two-point function reads [JT5], |Hj] 

S{t - t \x-x') = (<f>(x)cj>(x')) d = T J ^ e *^)-L cos^t -*'))• (31) 

For higher n-point functions one has 

(0 2 (xi) . . . 2 \x 2n )) d = 2 n S 2 (x 1 — x 2 ) ■ ■ ■ S 2 (x 2n -i - x 2n ) + permutations . (32) 

From Fig. 1 it is clear that the response function obtained after thermal averaging 
with Hq satisfies a Dyson-Schwinger equation (cf. Fig. 2), 

Di(x - x') = D R (x d 4 yd 4 y'D R (x - y)Tl(y - y') d D^{y' - x') . (33) 

The terms in the perturbative expansion for the self-energy II (x) are shown in Fig. 3 



+ 



Figure 2: Dyson-Schwinger equation for the retarded response function 

up to 0(X 2 ). The contribution O(X) (Fig. 3a) is linearly divergent and can be removed 
by a mass renormalization yielding the counter term [T(J 

The contribution 0(X 2 ) (Fig. 3b) is logarithmically divergent and can be rendered finite 
by a further mass renormalization 

Sml = -1a 2 /? / d 3 xS 3 (0,x) = --X 2 T 2 [ d$(0) — - — . (35) 
6 J 3 J ijj\ijj 2 ijjj, 



(a) (b) 

Figure 3: Perturbative expansion for the self-energy II. Dashed lines denote thermal two- 
point functions S. 

Contributions of higher order in A are all finite. 

We are interested in the damping rate, i.e., the imaginary part of the self-energy 
which is finite. Using Eqs. ([24]) and (|31|) one easily obtains for the self-energy to order 
A 2 , 

U(co, p) d = J dtd 3 x e ^e-^ U(t, x) d 

-A 2 T 2 Y m f d$(p) — , (36) 



Vl,V2,V3 



where rji = ±1, i = 1, . . . ,3. 

From Eq. (RBf) one obtains for the imaginary part of the self-energy (u> > 0), 



7r r I 

T(u>,p)d = -A 2 T 2 / d$(p) [uj^uj - a>i - a> 2 - toz) + ^<K^ + - a> 2 - io 3 )] . 

(37) 

We can now compare the classical damping rate T c i with the damping rate T of the 



full quantum theory as given by Eq. (14). At very high temperatures, i.e., T»m, the 



Bose-Einstein distribution function becomes 



f(u) ~ - » 1 . (3? 
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Inserting this approximation into Eq. (|T^) yields precisely Eq. (|37|). Hence, in the high- 
temperature limit, where number densities are very high, the damping rates are given 
by the classical theory to leading order. In particular the on-shell plasmon damping rate 
reads 

7 = h qu r{m, 0) = -i-/ gu A 2 T 2 , (39) 
2 153o7T 

where 

l gu = h(^T^ V2 (40) 

is the Compton wave length associated with the plasmon mass (D). Replacing l qu by l c \ 
yields the classical plasmon rate. This is the result of Aarts and Smit |TD[ . As our analysis 
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shows, this result extends to the full off-shell imaginary part of the self-energy in the 
high-temperature limit. 

Use of the Classical Theory 

So far we have performed the thermal average in the classical theory with respect to 
the free hamiltonian Ho, and we have calculated the self-energy II C ; to order A 2 . Higher 
orders in the expansion of Hr (cf. Fig. 1) and, furthermore, the thermal average with 
the full hamiltonian H rather than H yield corrections of higher order in A. 

However, if one is interested in the classical theory as a tool to evaluate the leading 
contribution of the quantum theory, the full result of the classical theory is irrelevant. 
Consider again the two-loop expression ( |Hp for T(uj,p). In the high-temperature limit 
terms proportional to the product of two Bose-Einstein distribution functions f{uji)f(ujj) 
dominate. This contribution is given by the classical theory. Contributions to T(u),p) 
proportional to a single distribution function /(o>j) are subdominant. Such terms do 
not occur in the classical theory. They lead to the following correction for the on-shell 
plasmon damping rate, 

8T(m, 0) = n\T f d$(0)— S(m + u t - u 2 - ou 3 ) , (41) 
J u 3 

which yields 

6 7 =^- 3 h\Hn(h\)T. (42) 

On dimensional grounds such a correction cannot appear in the classical theory. 

In summary, we have shown that in the high-temperature limit the imaginary part 
of the self-energy is determined by the classical theory to leading order in the coupling, 
whereas the dominant sub leading contribution is a genuine quantum effect. In order to 
obtain the leading order contribution to damping rates it is sufficient to perform the 
thermal average with respect to the initial conditions with the free hamiltonian. Such a 
simplification may be useful in numerical simulations for damping rates. 

We would like to thank D. Bodeker, A. Patkos and J. Polonyi for valuable discussions 
and comments. 



8 



References 

[1] V. A. Kuzmin, V. A. Rubakov and M. E. Shaposhnikov, Phys. Lett. B155 (1985) 
36 



[2] 
[3] 
[4] 
[5] 
[6] 
[7] 
[8 
[9 
[10 

[11 
[12 

[13 

[14 

[15 



D. Yu. Grigoriev and V. A. Rubakov, Nucl. Phys. B299 (1988) 67 

J. Ambjorn and A. Krasnitz, Phys. Lett. B362 (1995) 97 

G. D. Moore and N. Turok, preprint DAMTP 96-77, |hep-ph/ 9608350 



W. H. Tang and J. Smit, preprint ITFA-97-02, |hep-ph/9702"0T7| 
D. Bodeker, L. McLerran and A. Smilga, Phys. Rev. D52 (1995) 4675 

D. Bddeker, Nucl. Phys. B486 (1997) 500 
P. Arnold and L. McLerran, Phys. Rev. D36 (1987) 581 
P. Arnold, D. Son and L. G. Yaffe, preprint UW/PT-96-19, |hep-ph/ 9609481 

G. Aarts and J. Smit, Phys. Lett. B393 (1997) 395 
L. Dolan and R. Jackiw, Phys. Rev. D9 (1974) 3320 
R. R. Parwani, Phys. Rev. D45 (1992) 4695 

E. Wang and U. Heinz, Phys. Rev. D53 (1996) 899 

H. A. Weldon, Phys. Rev. D28 (1983) 2007 

G. Parisi, Statistical Field Theory (Addison- Wesley, 1988, New York), p. 309 



9 



